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The  paper  presents  a  framework  for  set-based  design  under 
uncertainty  and  demonstrates  its  viability  through  a  shape 
design  case  of  an  efficient  super-cavitating  hydrofoil  for 
ultra-high  speed  maritime  vehicles.  Redirecting  the  effort 
away  from  optimal  designs  to  those  that  safely  meet  relevant 
requirements  as  quantified  precisely  by  superquantile  mea¬ 
sures  of  risk  (s-risk),  the  paper  addresses  situations  faced  in 
practice  where  elimination  of  poor  designs  is  more  impor¬ 
tant  than  identifying  a  best  one.  The  complexity  of  the  design 
case,  involving  multiphysics  and  several  design  constraints, 
necessitates  the  use  of  surrogate  models,  which  are  built  us¬ 
ing  high-  and  low-fidelity  simulations  and  are  here  tuned  for 
the  first  time  to  reflect  s-risk.  Accounting  for  both  parameter 
uncertainty  as  well  as  errors  in  surrogate  models,  we  obtain 
robust  designs  that  satisfy  requirements  under  a  variety  of 
manufacturing  and  operating  conditions. 

1  Introduction 

It  is  well  known  that  an  optimized  design  tends  to  per¬ 
form  poorly  under  conditions  that  deviate  from  those  as¬ 
sumed  during  the  optimization.  In  this  paper,  we  propose  a 


two-pronged  approach  that  addresses  this  situation  and  leads 
to  designs  that  perform  robustly  under  a  variety  of  condi¬ 
tions.  The  first  prong  consists  of  redirecting  the  effort  away 
from  “optimal  designs”  and  instead  to  focus  on  identify¬ 
ing  sets  of  candidate  designs  that  meet  given  requirements. 
This  notion  of  set-based  design  was  pioneered  by  Toyota  and 
adopted  by  the  U.S.  Navy  [1],  It  responds  to  most  real-world 
design  processes  where  there  are  no  fixed  objective  and  con¬ 
straint  functions,  and  therefore  no  well-defined  optimal  de¬ 
sign  problem,  but  rather  an  evolving  collection  of  require¬ 
ments.  The  goal  is  then  to  find  designs  that  can  be  elimi¬ 
nated  from  further  consideration  and,  complementarily,  find 
the  candidate  designs  that  should  be  retained.  The  second 
prong  of  our  approach  is  to  rigorously  account  for  uncer¬ 
tainty  in  quantities  of  interest  (Qol)  and  their  models  in  a 
manner  that  captures  decision  makers’  preferences.  Here  we 
rely  on  superquantile  measures  of  risk  (s-risk)  [2 — 4]  to  as¬ 
sess  whether  Qols,  even  in  worst-case  scenarios,  meet  the 
specified  requirements. 

There  is  an  extensive  literature  on  design  under  uncer¬ 
tainty;  see  for  example  the  special  issue  [5]  and  the  recent 
articles  [4,  6-10],  Overviews  of  earlier  studies  are  found 
in  [11, 12],  It  is  beyond  the  scope  of  this  paper  to  provide 
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a  comprehensive  review  and  comparison.  However,  we  note 
that  much  of  the  literature  quantifies  performance  relative  to 
a  threshold  by  the  corresponding  failure  probability.  As  laid 
out  in  [13],  this  measure  is  fundamentally  problematic  as  it  is 
hard  to  estimate,  introduces  complexity  during  optimization, 
and  fails  to  account  for  the  severity  of  a  failure.  The  use  of 
s-risk  improves  on  these  aspects  as  described  below. 

For  the  first  time  we  consider  s-risk  in  a  real  design 
case  consisting  of  a  new  family  of  dual-operating  mode 
super-cavi taring  hydrofoils  [14]  devised  to  ensure  high  ef¬ 
ficiency  both  in  super-cavitating  and  in  fully-wet  conditions, 
unlike  the  conventional  geometries.  These  unconventional 
hydrofoils  are  designed  to  serve  as  lifting  surfaces  in  a  new 
ultra-high  speed  Hydrofoil  Small  Waterplane  Area  Twin  Hull 
(SWATH)  unmanned  craft  [15]  designed  to  operate  as  a  dis¬ 
placement  vessel  at  low  speeds  in  high  sea  states  and  to  reach 
a  top  speed  of  more  than  100  knots  in  milder  sea  states.  It 
is  evident  that  the  fully  wet  operation  of  the  unconventional 
hydrofoils  becomes  very  attractive  during  take-off  maneuver. 
In  this  design  case  and  in  fact  in  most  real-world  situations 
the  Qols  cannot  be  computed  exactly.  Numerical  simula¬ 
tions  allow  us  to  estimate  the  Qols  and  even  that  only  at  high 
computational  costs.  We  rely  on  surrogate  models,  trained 
by  output  of  the  simulations,  for  approximation  of  the  Qols 
across  the  space  of  design  variables.  Deviating  from  stan¬ 
dard  surrogate  models  (see  for  example  [16,  17]),  we  de¬ 
velop  for  the  first  time  a  risk-tuned  surrogate  model  that  is 
tailored  to  estimate  the  s-risk  of  Qols.  We  carry  this  out  in 
a  multi-fidelity  context  where  the  training  data  for  the  surro¬ 
gate  model  comes  from  a  low-fidelity  potential  flow  solver 
and  from  a  high-fidelity  viscous  fully  turbulent  model. 

The  contribution  of  the  paper  is  therefore  three-fold. 
First,  relying  on  s-risk,  we  develop  a  framework  for  set-based 
design  under  uncertainty  in  Section  2.  Second,  we  construct 
a  novel  surrogate  model  leveraging  multi-fidelity  simulations 
in  Section  3.  Third,  we  demonstrate  the  approach  through  a 
complex  design  case  in  Section  4,  with  results  in  Section  5. 

2  Set-Based  Design  under  Uncertainty 

In  this  section,  we  give  a  mathematical  formulation  of 
the  problem  of  set-based  design  under  uncertainty.  The  situ¬ 
ation  is  the  one  usually  encountered  in  practice:  the  process 
of  design  advances  through  several  stages.  At  an  early  stage 
there  are  many  designs  and  concepts  on  the  table,  as  well  as 
a  number  of  Qols  and  associated  requirements.  The  goal  at 
such  a  stage  is  not  to  identify  an  “optimal  design,”  but  rather 
a  set  of  candidate  designs  that  satisfies  the  requirements  and 
thus  are  worthwhile  to  consider  further.  Consequently,  the 
stage  eliminates  designs  deemed  deficient  according  to  the 
current  requirements.  The  set  of  candidate  designs  advances 
to  the  next  stage  to  face  a  revised  collection  of  Qols  and  re¬ 
quirements.  The  subset  of  designs  meeting  these  revised  re¬ 
quirements  becomes  the  new  set  of  candidate  designs  and  the 
process  is  repeated  at  successively  more  advanced  stages  of 
design.  We  here  provide  the  mathematical  foundations  for 
the  problem  of  identifying  the  set  of  candidate  design  at  a 
stage.  Our  formulation  applies  to  any  stage  despite  the  fact 


that  the  stages  might  involve  rather  different  Qols  and  re¬ 
quirements.  We  do  not  address  the  issues  of  selecting  Qols 
and  requirements  as  these  depend  on  the  design  case.  Section 
4  provides  details  in  the  context  of  shaping  a  hydrofoil. 

2.1  Accounting  for  Uncertain  Parameters 

Suppose  that  at  the  current  stage  of  the  design  process 
there  are  K  Qols.  The  numerical  value  of  each  Qol  is  pa¬ 
rameterized  by  a  vector  x  containing  design  variables  and 
a  vector  v  consisting  of  uncertain  parameters.  The  design 
variables  specify  a  design,  for  example  the  shape  of  a  hydro¬ 
foil  as  seen  below.  The  goal  is  to  find  designs  x  that  meet  the 
requirements  of  the  current  stage.  The  collection  of  such  x  is 
the  set  of  candidate  designs.  The  uncertain  parameters  rep¬ 
resent  inputs  that  are  beyond  the  control  of  the  designer  and 
that  can  affect  the  performance  of  the  design  such  as  envi¬ 
ronmental  conditions,  manufacturing  errors,  and  variability 
in  strength  of  materials.  Their  values  are  not  known  at  the 
current  stage.  We  let  the  functions  g*(x,v)  give  the  (true) 
numerical  value  of  the  kth  Qol  under  design  x  and  parame¬ 
ter  v,  k  =  1,2, ..., K.  Of  course,  in  most  practical  situations 
these  functions  are  not  available  in  explicit  form,  which  we 
address  below.  Without  loss  of  generality,  we  assume  that 
low  values  of  the  Qols  are  preferred  to  higher  values  and  that 
the  requirement  associated  with  the  kth  Qol  is  that  g/i(x,\) 
should  not  exceed  a  number  ;y.  If  a  practical  situation  de¬ 
mands  high  values,  we  simply  replace  g*(x,v)  by  —g/<(x,\) 
and  r/.  by  — r*.  If  both  too  high  and  too  low  values  of  the 
kth  Qol  need  to  be  avoided,  then  in  addition  to  ,^  (x.  v)  we 
also  consider  —gt i(x,v).  The  former  needs  not  to  exceed  an 
upper-bounding  requirement  and  the  latter  not  to  exceed  the 
negative  of  a  lower-bounding  requirement. 

A  requirement  that  a  design  x  should  have  g*(x,v)  not 
exceeding  r k  is  fundamentally  ill-defined  because  the  values 
of  the  uncertain  parameters  v  are  not  known.  We  assume  that 
the  uncertainty  in  v  can  be  quantified  by  a  known  probabil¬ 
ity  distribution.  Typically,  such  distributions  are  estimated 
based  on  historical  data  about  the  parameters,  but  we  omit  a 
discussion  of  this  subject  here;  see  for  example  [18].  When 
viewing  the  uncertain  parameters  as  a  random  vector  with  a 
given  distribution,  we  denote  it  by  V.  A  realization  of  V  is 
still  denoted  by  v.  Now  gk(x.  V)  is  a  random  variable  for 
each  x,  which  allows  us  to  bring  in  s-risk. 

S-risk  [13]  reduces  a  random  variable  to  a  number, 
which  somehow  represents  the  random  variable  and  can  be 
used  in  comparison  with  requirements.  It  was  originally  pro¬ 
posed  for  finance  applications  under  the  name  conditional 
value-at-risk  [2,3].  Specifically,  for  risk-averseness  param¬ 
eter  a  £  [0, 1],  the  s-risk  of  g^(x.V)  at  level  a,  denoted  by 
Ra(gk(x,V)),  is  understood  as  the 

average  value  of  gk(x,Y)  over  the  worst  (1) 

(1  —  a)100%  outcomes  of  g*(x,  V). 

This  expression  can  be  taken  as  the  definition  if  gk(x,Y) 
has  a  continuous  distribution  function.  If  the  risk-averseness 


parameter  a  =  0,  then  Ra(gk(x,Y))  is  simply  the  mean  of 
gk(x.Y).  If  a  =  1,  then  Ra(gk(x,Y))  is  the  largest  possi¬ 
ble  value  of  gk{x,  V).  Of  course,  a  value  of  a  between  these 
two  extremes  provides  a  middle  ground.  If  the  random  vari¬ 
able  gk(x,Y)  is  normal  with  mean  p  and  standard  deviation 
a,  then  Ra(gk(x,Y))  =  ju  +  G<p(0~1(a))/(l  —  a),  where  cp 
is  the  standard  normal  probability  density  function  and  <f>  is 
the  standard  normal  cumulative  distribution  function.  We  re¬ 
fer  to  [4]  for  general  methods  for  computing  s-risk.  Regard¬ 
less  of  the  situation,  s-risk  is  fundamentally  appealing  as  it 
has  a  simple  interpretation,  accounts  for  the  tails  of  proba¬ 
bility  distributions,  and  introduces  in  some  sense  no  signif¬ 
icant  complexity  compare  to  a  problem  without  uncertainty. 
Specifically,  if  g*(x,v)  is  convex  in  x,  then  Ra(gk(x,Y)) 
will  also  be  convex  in  x.  This  implies  that  finding  x  with 
^afe(xi^))  <  rk  is  in  some  sense  no  harder  than  finding  an 
x  with  gk(x,\)  <  rk  for  a  fixed  v.  In  contrast,  finding  an  x 
such  that  the  (failure)  probability  of  gk(x,  V)  exceeding  r 'k  is 
fundamentally  hard.  This  probability  also  fails  to  account  for 
the  magnitude  of  any  exceedance  of  rk\  see  [4, 13], 

We  are  now  in  a  position  to  formalize  the  requirement 
that  the  value  of  the  kth  Qol  should  not  exceed  rk.  Specifi¬ 
cally,  we  say  that  gk(x,  V)  is  safely  <  rk  if  Rafe(x,  V))  <  rk. 
In  view  of  Eqn.  (1),  it  is  clear  that  this  interpretation  of 
the  requirement  implies  that  even  on  average  over  a  set  of 
worst  outcomes,  the  value  of  the  Qol  will  not  exceed  rk.  The 
specific  definition  of  “worst  outcomes”  depends  on  a.  We 
can  also  interpret  the  requirement  probabilistically.  If  we  let 
the  event  that  gk(x,Y)  exceeds  rk  be  define  as  failure,  then 
^al?/t(x,  V))  <  rk  is  equivalent  to  having  the  buffered  fail¬ 
ure  probability  <  1  —  a  as  we  see  from  [13].  That  reference 
also  establishes  that  Ra(gi t(x,  V))  <  rk  implies  that  the  fail¬ 
ure  probability  is  <  1  —  a.  Letting  ak  £  [0,1]  be  the  risk- 
averseness  parameter  associated  with  the  kth  Qol,  the  above 
discussion  leads  to  the  following  requirements  on  the  design 
x: 

Rak fe(x,V))  <  rk  for  all  k  =  1,2 

This  formulation  accounts  for  the  uncertainty  in  the  parame¬ 
ters  V.  We  build  on  this  as  we  proceed  next  to  also  account 
for  uncertainty  in  surrogate  models. 

2.2  Accounting  for  Uncertainty  in  Surrogate  Models 

When  the  probability  distribution  of  gk(x,Y)  is  known, 
it  is  straightforward  to  compute  Rak(gk(x,Y))',  see  [4].  How¬ 
ever,  this  is  rarely  the  case  as  real-world  situations  involve 
most  often  complex  functions  gk(x,  v).  In  particular,  the 
probability  distribution  of  gk(x,  V)  is  usually  unknown  even 
if  that  of  V  is  given.  To  overcome  this  difficulty  we  utilize 
surrogate  models. 

S-risk  of  gi i(x,  V)  is  a  function  of  the  design  x,  which 
we  denote  by  sa(x),  i.e.,  sk(x)  =  Rak(gi i(x,V)).  When  this 
function  is  inaccessible,  as  is  usually  the  case,  we  develop  a 
surrogate  model  of  sk (x)  denoted  by  Sk(x).  Thus,  in  some 
sense  sk(x)  ss  Sk{x).  We  permit  that  Sk  is  a  random  field 


over  the  space  of  design  variables,  with  the  randomness  aris¬ 
ing  from  the  uncertainty  associated  with  the  accuracy  of  the 
surrogate  model.  Examples  of  surrogate  models  in  the  form 
of  random  fields  of  this  kind  are  those  based  on  kriging  and 
co-kriging;  see  for  example  [16, 17].  In  the  next  section  we 
provide  another  class  of  surrogate  models  that  are  also  of  this 
form. 

The  original  requirement  sk(x)  <  rk,  which  is  usually 
impossible  to  evaluate,  is  now  replaced  by  one  involving  the 
surrogate  model.  In  a  situation  analogous  to  that  above  we 
find  that  requiring  Sk(x)  <  rk  is  ill-posed  as  Sk(x)  is  a  random 
variable  for  every  x.  Thus,  again,  we  turn  to  s-risk.  For  risk- 
averseness  parameters  Pa  £  [0,1],  k  -  1,2,..., A,  we  adopt 
the  requirements 

R$k{Sk{x))  <  rk  for  all  k=  1,2,..., A. 

This  is  a  well-posed  requirement  with  the  interpretation  that 

Rpk  (Sk(x))  =  the  average  value  of  Sk (x)  over  the  worst 
(1  —  PQ 100%  outcomes  of  Sk(x). 

We  stress  that  here  the  average  is  over  possible  values  of  an 
uncertain  surrogate  model.  These  values  are  unrelated  to  the 
uncertain  parameters  V.  Since  the  surrogate  model  is  im¬ 
perfect,  Sk  (x)  is  random  and  we  account  for  this  fact  when 
making  the  comparison  with  rk. 

We  are  then  in  a  position  to  define  the  set  of  candidate 
designs  at  the  current  stage: 

designs  x  with  Rpk(Sk(x))  <  rk  for  all  k  =  1,2, ..., K.  (2) 

Designs  that  meet  these  requirements  are  safeguarded 
against  poor  performance,  accounting  for  both  parametric 
uncertainty  in  V  (using  risk-averseness  parameters  ak)  and 
surrogate  model  imprecision  (with  risk-averseness  parame¬ 
ters  Pa).  Designs  that  do  not  meet  these  requirements  are 
eliminated  from  the  design  process.  The  set  of  candidate 
designs  are  carried  forward  to  the  next  design  stage,  where 
they  will  faced  revised  requirements  for  possibly  new  Qols. 
The  process  is  then  repeated,  possibly  with  revised  risk- 
averseness  parameters. 


3  Risk-Tuned  Surrogate  Models 

The  previous  section  permits  any  surrogate  model  Sk(x). 
Here  for  the  first  time  we  develop  a  surrogate  model  that  is 
tuned  to  estimating  s-risk;  its  theoretical  foundation  is  given 
in  [19];  see  also  [20].  For  simplicity,  we  omit  the  subscript  k 
as  the  process  will  be  identical  for  each  Qol. 

For  a  Qol  given  by  g(x,v),  typically  from  a  high- 
fidelity  simulation,  we  leverage  a  low-fidelity  approxima¬ 
tion  /z(x,v).  There  is  no  difficulty  with  including  many 
levels  of  fidelity,  but  we  omit  the  details.  We  are  given 


a  data  set  {(g(,/z,-,x(-,  v,)  Ml  1  consisting  of  designs  x,-,  pa¬ 
rameter  values  v(,  and  the  corresponding  simulation  output 
gi  =  g(x/,Vi)  and  h,  =  /z(x,-,  v,).  We  adopt  a  model  such  as 
co  +  ci/z;  +C2hf  +  cTx,-  to  approximate  gi,  where  co,  ci,  and 
C2  are  scalar  coefficients  and  c  is  a  vector  of  coefficients  to 
be  determined  by  a  particular  type  of  error  minimization.  Of 
course,  models  with  more  or  less  terms  can  easily  be  con¬ 
sidered.  The  generation  of  the  data  set  requires  high-  and 
low-fidelity  simulations  of  the  relevant  Qol. 

The  essence  of  the  surrogate  model  is  the  method  by 
which  the  coefficients  in  the  model  are  fitted.  Specifically,  let 
a  €  (0, 1 )  be  the  risk-averseness  parameter  of  interest.  Then, 
the  proposed  method  stipulates  that  the  coefficients  should 
be  determined  through  a-quantile  regression,  i.e.,  find  coef¬ 
ficients  co,  ci,  C2,  and  c  that  minimize 

1  (  N  \  1  N 

(T^(g>ax{0' 

where  the  error  e,-  =  gi  —  (co  +  ci/z,-  +  C2hj  +  cTx,).  We 
discard  the  optimal  co,  but  name  the  other  optimized  co¬ 
efficients  ci,  C2,  and  c.  Finally,  we  set  co  =  the  s-risk  of 
{gi  —  ( c\hj  +  C2I1}  +  cTx,)}f=1.  Corollary  4.2  in  [19]  states 
that 

Ra{{gi}f=  1)  <f?a({co  +  Ci/z;  +  C2/z-  +CTX,}^=1). 

Consequently,  over  the  given  data  set  the  model  co  + 
ci/z(x,v)  +C2/z2(x,v)  +ctx  leads  to  an  upper  approximation 
of  the  s-risk  of  the  Qol.  This  motivates  the  use  of  the  model 
as  an  approximation  of  the  s-risk  of  the  Qol,  i.e., 

Ra{g{x,Y))  «/?a(c0  +  ci/z(x,V)+C2/z2(x,V)+cTx).  (3) 

There  are  two  reasons  for  why  this  approximation  might  fall 
short  of  an  actual  inequality.  First,  fitting  of  the  coefficients 
of  the  surrogate  model  might  be  carried  out  using  a  data  set 
of  size  N  and  not  the  true  distribution  of  the  random  vector 
V.  Second,  Eqn.  (3)  makes  comparisons  at  individual  de¬ 
signs  x,  which  might  generally  be  different  than  considering 
“averages”  across  a  variety  of  designs  as  relied  on  while  fit¬ 
ting  the  coefficients. 

We  note  that  it  is  central  that  the  risk-averseness  param¬ 
eter  a  used  in  the  regression  matches  that  of  interest  so  that 
the  coefficients  in  the  model  are  “tuned”  appropriately.  In  the 
case  of  a  =  0,  we  simply  determine  the  coefficients  by  means 
of  standard  least-squares  regression,  and  for  a  =  1  see  [19]. 

It  would  also  be  prudent  to  account  for  the  possible  inac¬ 
curacy  in  the  approximation  in  Eqn.  (3).  We  achieve  this  by 
replicating  the  above  estimation  on  different  data  sets  such 
that  not  only  one  but  J  sets  of  coefficients  are  determined. 
For  j  =  1  let 

S^(x)  =  Ra  (c^  +c[j\(x,V)  +  C2^h2(x,V)  +  (c«)Tx) 


be  the  resulting  models.  Finally,  we  define  the  surrogate 
model  as  S(x)  =  { S ^  (x)  H=| ,  which  can  be  viewed  as  a  ran¬ 
dom  field  over  the  space  of  design  variables,  with  J  possible 
and  equally  likely  realizations.  As  stipulated  in  Eqn.  (2),  a 
risk-averseness  parameter  p  £  [0, 1]  and  a  corresponding  s- 
risk  lead  to  the  condition  /?p(S(x))  <  r  for  designs  to  satisfy, 
where  r  is  a  requirement. 

In  summary,  we  see  that  through  Eqn.  (3),  the  vari¬ 
ous  (x)  approximate  the  s-risk  of  the  quantity  of  interest 
Ra(g(x.\)).  However,  to  account  for  the  fact  that  these  ap¬ 
proximations  are  uncertain,  we  consider  the  average  of  the 
worst-case  values  of  (x)  in  some  sense  that  depends  on 
p.  That  average  is  compared  with  the  given  requirement.  For 
example,  if  P  =  1,  then  S(x)  =  max7=i S^(x)  and  the 
worst-case  model  is  used  in  comparison  with  r.  For  lower  p, 
an  average  of  the  higher  (x)  would  be  used. 

4  Hydrofoil  Design 

As  a  practical  application  of  the  proposed  method,  we 
selected  the  design,  by  numerical  simulation,  of  a  new 
surface-piercing  super-cavitating  hydrofoil.  The  perfor¬ 
mances  achievable  with  a  CFD  based  optimization  offer  hy¬ 
drofoil  crafts  new  possibilities  of  reaching  higher  speeds  at 
sea  avoiding  the  hysteresis  phenomena  typically  affecting 
conventional  super-cavitating  profiles.  The  innovation  of 
this  hydrofoil  lies  in  the  shape  of  the  2D  sections,  which 
are  designed  to  ensure  particular  hydrodynamic  properties. 
Traditional  super-cavitating  hydrofoils  are  designed  to  oper¬ 
ate  at  very  high  speeds  only  where  flow  conditions  ensure 
the  presence  of  a  stable  vapor  cavity  enveloping  the  entire 
suction  surface  and  closing  many  chords  aft  a  blunt  trailing 
edge.  The  shape  of  the  unconventional  hydrofoil  presented 
in  [14,21]  features  a  pointed  leading  edge  ensuring  cavity 
detachment  at  the  operating  angle  of  attack  and  cavitation 
index  in  design  conditions  and  a  sharp  edge  on  the  pressure 
side  triggering  base  cavitation  when  working  at  zero  angle 
of  attack.  The  blunt  trailing  edge,  typical  of  conventional 
super-cavitating  profiles,  is  instead  tapered  in  a  tail  which 
is  designed  to  be  enclosed  in  the  supercavity  at  high  speed 
while  producing  a  good  pressure  recovery  and  higher  lift  at 
lower  speeds  when  cavitation  disappears.  This  guarantees 
the  performance  of  the  hydrofoil  also  in  sub-cavitating  con¬ 
ditions.  The  tail  is  functionally  separated  by  the  forward 
body  of  the  profile  though  a  sharp  corner  (face  cavitator)  on 
the  face  in  order  to  trigger  base  cavitation  at  intermediate 
speeds  or  at  high  speeds  and  lower  angles  of  attack.  Control 
points  of  composite  B-Spline  curves  provide  a  fully  paramet¬ 
ric  description  of  the  shape. 

A  benchmark  hydrofoil  has  been  designed  by  means  of 
a  differential  evolution  optimization  algorithm  (DE)  in  [22], 
where  CFD  simulations  were  performed  using  a  low-order 
Boundary  Element  Method  (BEM)  to  solve  for  the  steady 
potential  flow  with  cavitation.  Penalty  terms  were  added  to 
the  objective  function  when  design  variables  and/or  Qols  fell 
outside  of  given  bounds.  No  uncertainty  was  considered  in 
the  development  of  the  benchmark  design.  Figure  1  presents 
a  time-averaged  snapshot  of  the  flow  around  the  benchmark 


Fig.  1.  Time-averaged  vapor  content  contours  for  the  benchmark 
hydrofoil  resulting  from  an  optimization  without  accounting  for  uncer¬ 
tainty.  Flow  follows  positive  y  direction  (left  to  right).  Red  colors 
indicate  stable  cavity  enveloping  the  suction  surface  and  detaching 
from  the  face  cavitator  (sharp  edge  in  the  lower  hydrofoil  surface) 
with  Reynolds  number  4.58  •  107,  angle  of  attack  6  degrees,  and 
cavitation  index  0.05. 

design.  The  cavity  shape  can  be  inferred  from  the  color  scale 
indicating  vapor  content.  The  optimization  leads  to  a  thin 
cavity  enveloping  the  hydrofoil  surfaces  showing  the  draw¬ 
back  of  the  design:  it  is  highly  sensitive  to  changes  in  oper¬ 
ating  conditions.  In  this  case  the  existence  of  a  super-cavity 
relies  on  the  presence  of  a  vapor  regime  over  the  hydrofoil 
suction  surface.  It  is  evident  that  manufacturing  errors  or 
variable  conditions  (especially  in  terms  of  angle  of  attack) 
could  significantly  affect  the  flow  regime  and  hydrofoil  per¬ 
formances  for  a  design  with  a  thin  cavity.  This  recognition 
leads  us  to  considering  the  set-based  design  framework  un¬ 
der  uncertainty  described  above.  The  remainder  of  the  paper 
demonstrates  that  this  framework  produces  robust  designs 
that  perform  well  under  a  range  conditions. 

4.1  Description  and  Parameters 

We  provide  here  a  detailed  description  of  the  hydrofoil 
design  problem.  The  design  vector  x  consists  of  15  variables 
representing  coordinates  of  the  control  polygon  of  the  com¬ 
posite  B-spline  curves  describing  the  shape  of  a  hydrofoil. 
The  configuration  of  the  control  points  has  been  chosen  in 
such  a  way  that  all  desired  shape  variations  are  allowed  both 
on  the  suction  and  pressure  side.  Figure  2  gives  a  schematic 
representation  of  the  15  variables  together  with  a  resulting 
hydrofoil  shape.  The  pressure  side  is  described  by  five  con¬ 
trol  points,  one  of  which  allows  to  change  the  vertical  posi¬ 
tion  of  the  face  cavitator,  i.e.,  the  sharp  edge  at  the  beginning 
of  the  hydrofoil  tail.  Variations  of  the  suction  side  are  per¬ 
formed  by  five  control  points;  three  of  them  move  together 
vertically  following  the  coordinate  (xiQ.  The  leading  and 
trailing  edge  position  as  well  as  the  control  points  of  the  tail 
pressure  side  are  kept  fixed.  The  control  points  on  the  hydro¬ 
foil  face  regulate  lift  and  drag  performance  while  those  on 
the  back  suction  side  are  needed  to  control  the  cavity  thick¬ 
ness  and  maintain  a  sufficient  strength  (inertia  modulus). 

We  consider  uncertainty  related  to  manufacturing  errors. 
Specifically,  a  random  vector  V  =  (Vi , ...,  V15)  describes  rel¬ 
ative  displacement  values  to  the  control  point  coordinates. 


Table  1 .  Benchmark  design  and  bounds  of  design  variables. 


Variable 

Benchmark 

vmin 

Ai 

rmax 

Ai 

*1 

5.556- 

10^2 

5.556 -10^2 

1.667- 

10-1 

X2 

2.227  • 

10-1 

2.222- lO^1 

3.333- 

10-1 

X3 

3.900- 

lO^1 

3.889 -lO^1 

5.000- 

10-1 

X4 

5.576- 

10-1 

5.556 -lO^1 

6.667- 

10-1 

X5 

5.556- 

10^3 

0 

5.556- 

10-3 

X6 

1.111  • 

10^2 

1.111  - 10  3 

1.111- 

10-2 

XL 

3.456  • 

10^2 

1.333 -10^2 

3.756- 

10-2 

*8 

1.833- 

10^2 

5.556- 10^3 

2.778- 

10^2 

Xg 

5.556- 

10-3 

—  1.111  - 10  3 

5.556- 

10-3 

XlO 

1.979- 

10-1 

1.977  -lO^1 

2.417- 

10-1 

xn 

4.208  • 

10-1 

4.000-  lO^1 

4.889- 

10-1 

Xl2 

5.604- 

10-1 

5.500-  lO^1 

6.722- 

10-1 

Xi3 

3.833- 

10^2 

3.200- 10^2 

3.911- 

10^2 

*14 

6.711- 

10^2 

5.600- 10^2 

6.844- 

10^2 

*15 

8.267  • 

10^2 

7.044 -10^2 

8.267- 

10-2 

That  is,  the  actual  shape  manufactured  from  a  (nominal)  de¬ 
sign  x  =  (xi,  ...,^15)  is  assumed  to  be  x,-  =  x,-(l  +  Vj)  for  all 
i  =  1 , 15.  We  refer  to  the  vector  x  =  (xi,...,xis)  as  the 
(random)  shape  of  the  hydrofoil  under  design  x.  The  de¬ 
signer  selects  the  design  x  from  an  allowable  range  given  by 
the  bounds  x|llln  <  x,-  <  x',lax  for  all  i  =  1, ...,  15;  see  Table  1, 
which  also  gives  the  control  point  values  of  the  benchmark 
design  in  Figure  1.  The  values  of  the  design  variables  are 
relative  to  the  hydrofoil  chord. 

The  distribution  of  V  is  discrete  with  898  equally  likely 
realizations.  These  realizations  are  generated  by  sampling 
5000  15-dimensional  vectors  from  a  uniform  distribution 
on  [—0.06, 0.06] 15  and  then  eliminating  realizations  that  are 
deemed  unrealistic.  Thus,  “noise”  up  to  ±6%  is  added  to 
Xj.  The  framework  for  set-based  design  under  uncertainty  is 
indifferent  to  the  distribution  of  V.  Here,  we  select  a  dis¬ 
crete  distribution  to  avoid  errors  due  to  sampling  and  other 
approximations,  and  more  easily  highlight  the  features  of  the 
framework.  Using  chord  length  c  =  0.66m,  we  define  five 
(K  =  5)  Qols:  the  first  being  the  profile  inertia  w  that  is  re¬ 
quired  to  exceed  the  minimum  value  wmin  =  8.1  •  1 0  6  m3. 
The  second  Qol  is  the  profile  thickness  tj,  at  2%  of  the  chord 
that  needs  to  exceed  the  minimum  value  f™n/c  =  0.2%.  The 
third  Qol  defines  hydrofoil  performances  in  terms  of  drag- 
over-lift  ratio  Cd / 67,  that  needs  to  meet  a  requirement  of  0. 1 , 
where  Co  =  D/(0.5pcU2)  and  C/  =  L/(0.5pcU2)  normal¬ 
ize  drag  D  and  lift  L  under  operating  speed  U  =  61.667  m/s 
and  water  density  p  =  997.3  kg/m3.  The  fourth  Qol  is  lift 
Cl  having  a  requirement  of  0.2625  at  the  design  angle  of 
attack  CLdes  =  6  degrees  and  at  the  design  cavitation  num¬ 
ber  Odes  =  0.05.  The  fifth  Qol  is  that  the  negative  lift  —Cl 


Fig.  2.  Definition  of  design  variables,  horizontal  and  vertical  arrows 
indicate  the  allowed  movement  of  each  control  point,  and  correspond 
to  design  variables  x  =  (x\,  ...,xi5).  Some  additional  control  points 
are  fixed  and  illustrated  by  green  squares.  Three  control  points  on 
the  upper  suction  side  move  together  according  to  A15. 

must  meet  the  requirement  —0.23.  The  last  two  Qols  im¬ 
plements  the  condition  that  Cl  G  [0.23,0.2625].  The  risk- 
averseness  parameters  are  oti  =  0C2  =  0.7  (for  the  geometric 
properties),  (X3  =  0.9  (for  drag-over-lift),  and  C/4  =  (X5  =  0 
(for  lift).  Consequently,  we  seek  designs  that  on  average  in 
the  worst  30%  of  the  realizations  of  V  has  —w  <  —  wmm  and 
tj,  <  —t'pm.  Moreover,  on  average  in  the  worst  10%  of  the 
realizations  of  V  we  have  Co  /Cl  <0.1  and  the  average  of 
Cl  is  between  0.23  and  0.2625.  Although  suppressed  by  the 
notation,  w,tp,Co,  and  Cl  depend  on  the  shape  of  the  profile 
and  therefore  on  both  the  design  x  and  the  random  vector  V. 

The  first  two  Qols  will  not  need  surrogate  models  as 
there  are  explicit  formula  for  inertia  and  profile  thickness. 
For  drag-over-lift  and  lift  we  develop  surrogate  models  of  the 
kind  given  in  Section  3  using  P3  =  0.95  for  drag-over-lift  and 
p4  =  P5  =  0.  for  lift.  The  number  of  replications  of  surrogate 
models  is  set  to  J  =  9,  each  fitted  using  roughly  a  ninth  of  the 
data  as  described  below.  The  choices  of  risk-averseness  pa¬ 
rameter  p3  then  implies  that  the  value  of  the  worst  surrogate 
model  over  these  9  is  used  in  comparison  with  the  require¬ 
ment  of  0.1.  For  the  lift,  the  average  is  employed. 

4.2  Analysis  Models  and  Simulations 

We  have  to  our  disposal  a  high-fidelity  solver  for  an  un¬ 
steady  viscous  method  for  the  solution  of  Reynolds  Aver¬ 
aged  Navier-Stokes  equations  (URANSE)  as  well  as  a  low- 
fidelity  potential  flow  solver  that  both  return  for  a  given 
shape  x=  (jci, ...., *15)  the  lift  Cl  and  the  drag  Co  at  the  op¬ 
erating  conditions.  The  URANSE  and  potential  flow  output 
are  referred  to  as  high-  and  low-fidelity  output,  respectively. 
In  the  present  context,  each  high-fidelity  solution  takes  about 
6  hours  on  4  cores,  while  each  low  fidelity  run  requires  5  sec¬ 
onds  on  a  single  core.  Figure  1  displays  high-fidelity  output. 

The  potential  flow  solver  is  based  on  a  BEM  in  which 
a  Laplace  equation  is  solved  for  a  perturbation  potential  and 
a  Bernoulli  equation  is  solved  instead  of  continuity  and  mo¬ 
mentum  equations  for  a  steady  non-viscous  flow.  Friction 
forces  on  the  wetted  surface  of  the  hydrofoil  are  included  us¬ 
ing  an  empirical  local  friction  coefficient.  A  cavitation  model 
solving  for  the  cavity  thickness  on  the  hydrofoil  suction  sur¬ 
face  is  included  by  imposing  an  additional  set  of  kinematic 
boundary  conditions  that  are  satisfied  through  additional  sin¬ 
gularities  distributed  on  the  hydrofoil  suction  surface  and 


wake  panels.  This  simplified  model  is  able  to  predict  mid¬ 
chord  cavity  closure.  Since  face  cavitation  does  not  repre¬ 
sent  a  feasible  engineering  design,  bubble  dynamic  has  not 
been  modeled  for  cavities  developing  on  the  hydrofoil  pres¬ 
sure  surface;  see  [22]  for  a  detail  description  of  the  potential 
flow  solver  in  this  setting. 

The  fully  turbulent  solution  of  URANSE  is  achieved  for 
a  multi-phase  flow  consisting  of  a  fluid  mixture  of  liquid  and 
vapor.  The  problem  is  numerically  approached  using  a  fi¬ 
nite  volume  technique  with  collocated  arrangement  of  vari¬ 
ables  where  PDEs  are  solved  for  the  fluid  mixture  through 
a  scalar  quantity  y  indicating  the  relative  content  of  vapor 
with  respect  to  liquid  within  each  cell.  This  surface  captur¬ 
ing  method  developed  by  [23]  is  known  as  Volume  of  Fluid 
technique  and  it  is  particularly  suitable  for  the  solution  of 
complex  fluid  interfaces  as  the  one  characterizing  unstable 
cavitating  flows.  An  additional  transport  equation  for  y,  rep¬ 
resenting  the  vapor  content,  is  solved  together  with  RANS 
equations.  This  transport  equation  is  solved  once  the  source 
term  indicating  the  specific  net  mass  transferred  in  the  cavi¬ 
tation  process  is  known.  Condensation  and  vaporization  reg¬ 
ulating  the  above  mentioned  mass  transfer  are  found  through 
the  cavitation  model  formulated  by  [24],  In  this  simplified 
method  the  uncompressible  assumption  applied  to  the  liq¬ 
uid  phase  is  extended  to  the  vapor  phase,  hence  no  addi¬ 
tional  thermodynamic  equations  are  required  to  model  the 
phase  change  which  relies  only  on  the  local  pressure  within 
each  control  volume.  The  super-cavitating  flow  around  the 
hydrofoils  is  solved  using  a  hybrid  structured-unstructured 
grid  consisting  in  three  different  regions:  a  structured  grid 
is  used  close  to  the  hydrofoil  surface  to  accurately  predict 
the  strong  gradients  characterizing  boundary  layer  flows,  two 
nested  unstructured  grids  of  tetrahedral  elements  are  used  to 
increase  the  resolution  in  the  region  of  cavity  development 
(hydrofoil  wake)  while  keeping  a  coarse  discretization  in  the 
far  field  region;  see  [25]  for  details. 

For  the  shapes  resulting  from  the  benchmark  design  with 
the  898  realizations  of  V  we  run  the  URANSE  solver  and 
obtain  C/j  and  Cl  for  all  898  resulting  shapes.  The  same  in¬ 
put  to  the  potential  flow  method  results  in  some  estimates 
of  Co  and  Cl,  but  also  approximately  50%  computational 
failures,  due  to  the  predicted  presence  of  face  midcord  cav¬ 
itation  that  is  not  presently  resolved  by  the  method.  Conse¬ 
quently,  instead  of  using  the  actual  low-fidelity  simulations 
for  each  of  the  898  shapes,  we  compute  for  each  shape  a 
weighted  average  over  all  the  successful  runs,  with  weights 
determined  by  the  distance  from  the  current  shape  to  the 
shape  of  those  runs.  Throughout,  we  deal  exclusively  with 
this  modified  low-fidelity  data  set.  Figure  3  illustrates  the 
high-  and  low-fidelity  drag-over-lift  values  across  the  898 
shapes.  These  runs  generate  the  data  from  which  the  surro¬ 
gate  models  are  fitted  according  to  the  procedure  in  Section 
3.  Since  0C3  =  0.9,  the  surrogate  model  for  drag-over-lift  is 
achieved  by  carrying  out  quantile  regression  with  a  =  c*3. 
The  data  set  of  898  drag-over-lift  values  from  the  highl¬ 
and  low-fidelity  simulations  is  divided  into  7  =  9  roughly 
equally  large  groups,  on  which  quantile  regression  is  car- 
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Fig.  3.  Low-  and  high-fidelity  estimates  of  drag-over-lift.  Values  in 
meter. 
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Fig.  4.  Design  2  (blue  curve,  triangles  for  control  points)  compared 
with  benchmark  profile  (black  curve,  squares  for  control  points). 


ried  out  and  a  surrogate  model  Sp,L(x)  is  computed  for  each 
j  =  1 , . . . ,  9.  Since  we  would  like  to  account  for  the  variability 
in  the  models,  we  adopt  the  whole  ensemble  as  the  surrogate 
model  for  drag-over-lift,  i.e.,  SDjL{x)  =  {5'^i(x)}®=1  and 

%(Sd/l(x))  =  max;=l  .  .  9.S'^L(x)  because  P3  =  0.95. 

The  surrogate  model  for  lift  is  constructed  similarly. 
Since  0C4  =  as  =  0,  the  surrogate  model  for  lift  is  achieved 
by  carrying  out  least-squares  regression.  We  rely  on  the  same 
data  set  as  above,  which  again  is  divided  into  J  =  9  groups 
that  result  in  surrogate  models  S^\x),j  =  1,...,9,  and  the 
overall  model  Sl(x)  =  {S^(x)}?=1.  Since  p4  =  0,  this  re¬ 
sults  in  /?p4(Sz,(x))  =  £y=i  sj^(x)/9.  For  the  negative  lift, 
we  simply  reverse  the  sign. 


Fig.  5.  Design  4  (red  curve,  circles  for  control  points)  compared  with 
benchmark  shape  (black  curve,  squares  for  control  points).  Values  in 
meters. 


Fig.  6.  Time-averaged  vapor  content  (left  column)  and  pressure  co¬ 
efficients  (right  column)  contours.  From  top  to  bottom:  benchmark 
design,  Design  2,  and  Design  4.  Flow  conditions  as  in  Figure  1. 


5  Results  and  Candidate  Designs 

We  are  now  in  a  position  to  generate  sets  of  candidate 
designs.  For  the  five  Qols,  we  have  the  conditions 


Rai  (-w)  <  -H’mm 
Ra2(-t2P)  <  -tfn 
Rh(SD/L(x))<  0.1 
%(Sl(x))<  0.2625 
Rfe(-St(x))<-  0.23 


(sufficiently  large  inertia) 
(sufficiently  large  thickness) 
(sufficiently  low  drag-over-lift) 
(sufficiently  low  lift) 
(sufficiently  high  lift) 


defining  the  candidate  designs,  where  the  first  two  re¬ 
quirements  avoid  surrogate  models  and  thus  and  /x’p, 
become  immaterial.  (Since  Rai{—w)  is  deterministic, 
( Rai  (— w))  =  Rai  (— w)  and  similarly  for  the  second  re¬ 
quirement.)  A  design  x  that  satisfies  these  inequalities  is 
guaranteed  to  meet  the  specified  requirements  in  a  mathe¬ 
matically  precise  sense  accounting  for  both  parametric  un- 


certainty  in  the  manufacturing  process  of  the  hydrofoil  and 
uncertainty  in  surrogate  models. 

The  above  five  conditions  are  trivial  to  evaluate  for  each 
design  x  as  they  involve  only  explicitly  or  easy-to-compute 
quantities.  No  high-fidelity  simulations  are  required.  Low- 
fidelity  output  is  needed  to  evaluate  the  surrogate  models  due 
to  the  presence  of  h  in  Eqn.  (3),  but  they  are  quickly  ob¬ 
tained.  It  is  therefore  easy  to  check  whether  a  given  design 
satisfies  these  conditions  and  also  to  generate  a  collection 
of  candidate  designs  satisfying  the  conditions.  For  the  sake 
of  demonstration,  we  simply  randomly  generate  500  designs 
within  the  bounds  of  Table  1  and  for  each  check  the  five  con¬ 
ditions.  This  results  in  four  candidate  designs  that  satisfy 
the  conditions.  (Many  more  could  have  been  obtained  by 
further  sampling  and  other  methods.)  In  the  following,  we 
concentrate  on  two  of  them,  labelled  Design  2,  with  (in  me¬ 
ters)  x  =  (0.0834,  0.2629,  0.4258,  0.5611,  0.0040,  0.0099, 
0.0147,  0.0094,  0.0027,  0.1887,  0.3929,  0.5602,  0.0300, 
0.0602,  0.0742),  and  Design  4,  with  (in  meters)  x  =  (0.1069, 
0.2724,  0.3749,  0.5962,  0.0009,  0.0072,  0.0283,  0.0084, 
0.0045,  0.1877,  0.4047,  0.5058,  0.0297,  0.0578,  0.0742); 
see  Figures  4  and  5.  Designs  2  and  4  are  quite  different,  with 
the  two  other  candidate  designs  resembling  these. 

Table  2  summarizes  the  performance  of  Designs  2  and 
4.  Columns  2  and  3  give  the  drag-over-lift  for  Designs  2  and 
4,  respectively.  Rows  3-11  give  the  estimates  provided  by 
the  surrogate  models  S^jL(x),  j  =  1,  ...,9,  resulting  in  the  as¬ 
sessment  Rp3  (So/l(x))  given  in  the  second-to-last  row.  The 
latter  quantity  was  compared  with  the  requirement  0.1  to  de¬ 
termine  whether  the  designs  satisfied  the  requirement.  For 
the  sake  of  validation,  we  compute  the  actual  s-risk  of  drag- 
over-lift  using  898  high-fidelity  simulations  and  obtain  the 
results  in  the  last  row.  We  observe  that  although  some  indi¬ 
vidual  surrogate  models  (for  example  S ^  (x))  underestimate 
drag-over-lift,  Rp3(SD/L(x))  provides  conservative  estimates 
of  the  actual  drag-over-lift.  This  highlights  the  importance  of 
considering  the  uncertainty  in  surrogate  models.  Columns  4 
and  5  give  similar  results  for  lift;  again  surrogate  models  are 
reasonably  accurate.  In  comparison,  the  benchmark  design 
has  actual  s-risk  for  drag-over-lift  of  0.0880  and  for  lift  of 
0.2675,  the  latter  being  excessively  high. 

The  s-risk  of  the  inertia  for  the  benchmark.  Design  2, 
and  Design  4  are  7.400-  10^6,  8.738  •  10~6,  and  8.48489  • 
10  6,  respectively.  We  see  that  when  accounting  for  man¬ 
ufacturing  errors,  the  benchmark  design  fails  to  meet  the 
requirement  of  wmm  =  8.1  •  1 0  6 .  The  s-risk  for  the  cor¬ 
responding  profile  thicknesses  are  0.00118,  0.00135,  and 
0.00185.  Again,  the  benchmark  design  fails  to  meet  the  re¬ 
quirement  of  f™ln  =  0.00132. 

Figures  4  and  5  shows  that  the  suction  side  control  points 
of  Designs  2  and  4  have  moved  to  achieve  a  thicker  cavity  by 
lowering  the  back  profile  surface.  The  minimum  hydrofoil 
thickness  at  2%  of  the  cord,  required  to  avoid  fluttering,  has 
been  recovered  shifting  the  pressure  side  control  points  close 
to  the  leading  edge  of  the  hydrofoil.  The  control  points  of 
the  pressure  side  have  been  changed  in  different  ways  to  en¬ 
sure  the  overall  hydrodynamic  performance:  Design  2  has 


Table  2.  Performance  of  Designs  2  and  4. 


Drag/Lift 

Lift 

Model 

Design  2 

Design  4 

Design  2 

Design  4 

S(1)(x) 

0.0959 

0.0971 

0.2354 

0.2367 

S(2)(x) 

0.0937 

0.0949 

0.2346 

0.2368 

S(3)(x) 

0.0925 

0.0930 

0.2346 

0.2397 

S(4)(x) 

0.0958 

0.0988 

0.2341 

0.2355 

S(5)(x) 

0.0954 

0.0982 

0.2352 

0.2368 

S(6)(x) 

0.0954 

0.0990 

0.2323 

0.2337 

S<7)(x) 

0.0962 

0.0989 

0.2345 

0.2381 

S(8)(x) 

0.0932 

0.0951 

0.2345 

0.2381 

S<9)(x) 

0.0956 

0.0970 

0.2363 

0.2390 

*p(S(x)) 

0.0962 

0.0990 

0.2346 

0.2371 

Actual 

0.0956 

0.0971 

0.2396 

0.2465 

a  single  curvature  shape  with  a  maximum  at  about  x  =  0.42 
m,  while  Design  4  shows  general  higher  curvatures  with  a 
marked  slope  change  at  about  x  =  0.28  m.  Design  4  has 
pressure  side  shape  similar  to  that  of  the  benchmark  design, 
but  the  maximum  curvature  of  the  concave  part  is  shifted  to¬ 
wards  the  trailing  edge.  Design  2  has  a  single  curvature  face 
and  Design  4  a  double  curvature  face.  Design  2  is  similar  to 
the  Johnson  two-terms  profile  [26],  also  known  as  the  Tulin- 
Burkat  hydrofoil  [27];  the  single  curvature  pressure  surface 
leads  to  lower  lift  force  generation,  which  for  this  particular 
shape  has  been  compensated  through  the  vertical  displace¬ 
ment  of  the  face  cavitator.  Moving  the  sharp  edge  on  the 
pressure  side  has  the  direct  effect  of  changing  the  cavity  de¬ 
tachment  position  and  virtually  increasing  the  operating  an¬ 
gle  of  attack.  Design  4  is  characterized  by  a  double  curva¬ 
ture  resembling  a  Johnson  five-terms  profile  and  leading  to 
lower  drag-over-lift,  i.e.,  improved  efficiency.  As  shown  in 
Figure  5,  the  vertical  position  of  the  face  cavitator  does  not 
move  compared  with  the  benchmark  shape.  Hence,  the  in¬ 
crease  in  lift  generation  is  mainly  due  to  the  double  curvature 
leading  to  a  double  pressure  peak  as  evidenced  in  Figure  6, 
which  shows  pressure  as  well  as  vapor  content  for  all  there 
designs.  Design  4  experiences  a  second  pressure  peak  up¬ 
stream  the  face  cavitator,  locally  behaving  as  a  stagnation 
point.  The  presence  of  bump  on  the  face  in  between  the  two 
pressure  peaks  generates,  as  seen  in  Figure  6,  a  lower  pres¬ 
sure  region  within  the  maximum  curvature  points.  Similarly 
to  the  situation  for  five-term  profiles,  this  region  could  be 
of  interest  in  certain  flow  conditions.  Both  Designs  2  and  4 
present  a  reduced  suction  surface  curvature  with  respect  to 
the  benchmark  design.  This  feature  directly  affects  cavity 
thickness  close  to  the  leading  edge  as  seen  in  Figure  6  and, 
hence,  ensures  more  robust  cavitating  regimes  than  for  the 
benchmark  profile.  The  set-based  design  framework  under 
uncertainty  thus  generated  designs  that  have  slightly  lower 


efficiency  due  to  a  larger  cavity  thickness,  but  that  are  more 
robust  to  changes  in  geometry  due  to  manufacturing  errors. 
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